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SUMMARY 

A complete formulation and numerical results for the collisionless space-charge 
problem of a cylindrical diode with thermionically emitted electrons are presented. This 
formulation allows calculation of the current-voltage characteristics from the Schottky 
retarding region to saturation. In addition, the onset of the Schottky retarding region is 
unambiguously defined. 


INTRODUCTION 

The study of collisionless electron flow in cylindrical geometries has a long history; 
a review of these studies is to be found in the article by Ivey (ref. 1). Except for the 
studies of Wheatcroft (1940, ref. 2) and Schottky (1914, ref. 3), the early studies in- 
volved broad idealized assumptions such as monoenergetic electron emission and virtual 
cathodes at the potential minimum. 

Schottky (ref. 3) derived an expression for the current -voltage characteristic in the 
retarding potential range, that is, for cases where the potential is monotonically decreas- 
ing. He assumed, however, that his relation was valid for all monotonically retarding 
potentials. It will be shown in this report that there exists, in general, a range of col- 
lector potentials associated with monotonically retarding potentials for which Schottky's 
relation is not valid. In his defense, however, let it be noted that this potential range 
becomes negligible for large ratios of collector radius to emitter radius (e.g. , filament 
emitters). 

Wheatcroft (ref. 2) succeeded in obtaining numerical solutions to the mathematical 
equations over a limited range of collector potentials. These solutions were obtained by 
laborious numerical calculations before the advent of the high speed electronic computer. 
He was forced, therefore, to employ some approximate methods of solution which unfor- 
tunately hid some of the physical aspects in the range he treated. Because of this limi- 
tation, numerical solutions for the cylindrical diode from the Schottky retarding regime 



to saturation have never been obtained. This report is an attempt to fill the void. We are 
more fortunate than Dr. Wheatcroft in the availability of high speed electronic computers. 

Because the wide spread availability of computers today provides the means of calcu- 
lating diode characteristics for all collector potentials of interest, it will not be necessary, 
as in previous studies (refs. 2 and 4), to concern ourselves with the position and magnitude 
of the potential minimum in the space -charge -limited cases. Meaningful to the experi- 
mentalist are the parameters he can control and the variables which he can observe. 

Among the former are emitter radius, temperature, work function, collector radius, and 
collector potential. The easiest (and often times the only) variable he can measure is the 
current to the collector. It is therefore somewhat artificial to present theoretical results 
in terms of the position and magnitude of the potential minimum. For given physical 
diode characteristics, it is, furthermore, easy to investigate the theoretical variation of 
the potential minimum by numerical solutions. 


ANALYSIS 

An investigation of the characteristics of the collisionless cylindrical diode from the 
Schottky retarding region (hereinafter SRR) to saturation is presented. SRR is defined 
as that region in which the Schottky relation (ref. 3) 

— = ae^erf erfc(lA^) (1) 

*o a 


is valid. In equation (1), 



( 2 ) 


( 3 ) 


r and r are the emitter and collector radii, respectively, V is the potential energy 
(not just potential) of an electron in the electrostatic field of the diode, T is the emitter 
temperature, k is the Boltzmann constant, and q? is the dimensionless electrostatic 
potential. (All symbols are defined in appendix A. ) 
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Figure 1. - Space-charge potentials. 


Physically, SRR denotes a region wherein the current to the collector is a function of 
the collector potential only and not of the functional form of the potential in the inter - 
electrode space. In the planar case (ref. 5), the onset of SRR corresponds to the onset 
of a monotonically retarding potential. In general, this correspondence is not true in the 
cylindrical case; that is, there may be monotonically retarding potentials for which the 
electron current to the collector is not a function of the collector potential alone. 

Characteristic potentials (and their corresponding collector potential) for the cylin- 
drical diode are shown in figure 1. Here the dimensionless radius r is defined as 
f/r Q . In general, the potential signaling the onset of SRR lies below that of the first 
monotonically decreasing potential, as shown by curves c and d, respectively. Curve b 
depicts a typical space -charge potential and curve a signals the onset of saturation cur- 
rent to the collector. Curve a is identified by zero potential slope at the emitter, while 
curve c is identified by zero potential slope at the collector. 

The problem is to find self-consistent solutions to Poisson's equation which, for the 
cylindrical diode, can be expressed in the form (derived in appendix B) 

<?"(«) = B Q e 2 ^,cp) (B14) 


where 


£ = In r 

the space -charge parameter B Q (mks units) is given by 


(B13) 
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(B12) 




2 2 
e n o r o 

e o kT 


and n (%,<p) represents the dimensionless electron density. 

With the assumption that the emitted electrons have a half-Maxwellian velocity dis- 
tribution and encounter no collisions, formal expressions for the electron density 
n(£, cp) as a function of radius and potential are derived in appendixes B to D. These 
expressions for the density become quite complex because of the exchange between radial 
angular momentum in cylindrical geometries. This exchange can be clarified by a 
description of the possible electron trajectories in terms of an equivalent one -dimensional 
potential energy. 


Equivalent One-Dimensional Potential Energies 


The sum of the radial and tangential energy components of an electron at radius 
may be expressed by 


- mu 2 + V(r) + 
2 



( 4 ) 


where L is the angular momentum. (Here dimensional variables are used for heuristic 
reasons. ) Since this sum is a constant of the motion, 


1 mu 2 + V(r) + 

2 




+ 



( 5 ) 


where, for convenience, V{r Q ) is taken to be 0. A new constant of motion E will now be 

2 a o 

defined by adding -L /2mr‘ to both sides of equation (5): 


E s 


1 

2 


‘2 

mu 


+ V(r) + 



2m \ f2 



( 6 ) 


Note that E is equivalent to the radial kinetic energy of emission. Now, as in refer- 

r*-> 

ence 6, an equivalent one -dimensional potential energy (hereinafter EODP) V can be 
defined as 
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( 7 ) 


V(r) = V(r) + II. (± - Jl\ = v(f) + V T (?) 

2 m V 2 v 2 J 

where the force on the equivalent particle is given by -dV/dr. The ’’angular potential 
energy” V-^(r) represents the contribution of angular motion to V(r). The potential 
energies V(r), V L (r), and V(r) are shown in figure 2 for a typical space-charge con- 



Figure 2. - Equivalent one-dimensional potential for cylindrical diode. 


dition. Here, V(r) has the form of a potential energy barrier in front of the emitter, 
while V-j j (r) has the form of a potential energy sink. Their sum, V(r), also has the 
form of a potential barrier in front of the emitter, but its maximum is much less than 
that corresponding to V(r). By referring again to V(r), it can be observed that only 
electrons with E > Eg can ’’escape” the emitter and reach the collector. This is 
in contrast to the planar geometry where E must be greater than Eg for the electron 
to reach the collector. An electron with E = Ej < Eg will reach a radius r j before 
returning to the emitter. 

Qualitative characteristics of the different EODP appearing in the analysis of the 
cylindrical diode space-charge problem are shown in figures 3 and 4. In these figures, 
the potential energy V(r) is a characteristic of the diode operating condition, while 
Vj^r) is the angular potential energy of a particular electron. In other words, the whole 
spectrum of angular potential energies is present at any given time since 
V^(r) oc l = TqVq, where v Q is the tangential electron emission velocity. For obvious 
reasons of clarity, only one V^r) is depicted in each case. Since the different EODP 
depend on the functional form of V(r), not all the types shown in figures 3 and 4 can exist 
at one time. 
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Emitter Collector 
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Emitter 


Collector 

Emitter 


Collector 


(0 



(d) 



Figure 3. - Types of equivalent one-dimensional potentials for cylindrical diode with potential maximum (r m is 
position at potential maximum). 


In figure 3(a) only those electrons for which E is greater than the maximum V 
(hereinafter Iwill reach the collector. The EODP of figure 3(b) allows all elec- 

trons to reach the collector regardless of E (even electrons for which E = 1/2 mu^ = 0). 
Electrons for which the EODP is as that depicted in figure 3(c) will reach r j regardless 
of E, but will reach the collector only for E > V . In figure 3(d), two maximums oc- 
cur in V. Electrons for which E is greater than the first V may still be returned 
to the emitter if E is not simultaneously greater than the second V max . 

Three other EODP are shown in figure 4 for the case of a monotonically retarding 
potential energy. The interpretations of the three different EODP are the same as in 
figure 3. Figure 4 is important because it clarifies the definition of the Schottky retard- 
ing region (SRR). As previously defined, SRR occurs when the electron current to the 
collector is a function of the collector potential only. It has been observed that an elec- 
tron cannot reach the collector unless it has an E > Hence, the SRR is defined 

by those potential energies V such that for all V T , V mov occurs at the collector. 

l max 
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(a) 



(b) 



(c) 

Figure 4. - Types of equivalent one-dimensional 
potentials for cylindrical diode with monotonic 
potential. 


Figure 4(a) shows that this condition is not always satisfied by a monotonically retarding 
potential energy. 

As shown in appendix B, Wheatcroft (ref. 2) implicitly recognized only that condition 
depicted in figure 3(a). A complete formulation of the collisionless electron space-charge 
problem embodying all the trajectories is given in appendixes B to D. 
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DISCUSSION 


In the absence of space charge, the current-voltage curves would coincide with the 
Schottky curves (defined by eq. (1)) in the retarding region (fig. 5). For r c /r = 1 
(In r c /r Q = 0), the Schottky curve is simply the Maxwell-Boltzmann line associated with 
the planar diode. The current-voltage curves with space charge present are shown in 
figure 6 for various values of In r c /r Q and B Q (defined in eq. (B12)). 

For given In r c /r Q , increasing B Q results in a larger saturation potential and a 
decrease in I/I Q for a fixed collector potential. A similar observation is valid if B Q 
is given and In r c /r Q is allowed to increase. 



Figure 5. - Schottky curves. 
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Figure 6. - Current-voltage curves. 
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(b) In r c /r 0 = 1. 
Figure 6. - Continued. 
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Dimensionless collector potential, <pi r c ) 





Space-charge 

parameter, 



Dimensionless collector potential, p(r c ) 




Space-charge 

parameter, 
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Dimensionless collector potential, <p\r c ) 




As mentioned in the INTRODUCTION, the onset of the SRR does not coincide, in 
general, with the onset of monotonic retarding potentials. Figure 7 demonstrates this 
fact quantitatively. The solid lines in this figure are simply the Schottky curves of fig- 
ure 5 plotted against a different abscissa to effect a better separation. The dashed lines 
indicate the locus of points where the potential first becomes monotonic for the indicated 
values of B Q . The dot-dash lines represent the locus of points signaling the onset of the 
SRR for the indicated values of B Q . Here it is also apparent that the two lines coalesce 
into one as r /r gets very large. 

u U 



Figure 7. - Limits on Schottky curve. 
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Comparisons of the theoretical results with some experimental results are given 
in figure 8. (The abscissas have been corrected for the experimentally determined con- 
tact potential A cp . ) The experimental results in figure 8 were reported by Schottky 
(ref. 3) in 1914. The theoretical curves of figure 8 were obtained by employing the ex- 
perimentally determined values of I Q , T, and In r c /r to compute the value of B Q 
and kT/e. 

There is, of course, no means at this late date of interrogating the experimentalist 



(a) Schottky (ref. 3). Run 1; collector radius, 56. 



(b) Schottky (ref. 3). Run 2; collector radius, 56. 
Figure 8. - Comparison with experiment 
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(c) Schottky (ref. 3). Run 8; collector radius, 30. 

Figure 8. - Concluded. 

to try to discover reasons for the discrepancies between theory and experiment. It is 

to be noted that Schottky also observed marked discrepancies between observed results 

which he could not satisfactorily explain. Of importance here, however, is that the 

entire theoretical current-voltage characteristic can be generated and compared with 

experiment. After perusing the theoretical current-voltage curves of figure 6, one 

might suspect that something had affected the experimental determination of B Q in 

figure 8(b); while in figure 8(c) either r /r or a combination of work function and B 

CO o 

might be suspect. 

Further experimentation might also better define limitations of the model. For 
instance, just when it is no longer valid to neglect short range electron-electron inter- 
actions is still a moot point. 


CONCLUSIONS 

A complete formulation of the collisionless cylindrical diode with electron emission 
has been presented. The numerical solutions describe the physical characteristics of 
these diodes. It is now possible, by comparing theory with experiment over the whole 
current-voltage curve, to obtain a more detailed understanding of the cylindrical elec- 
tron diode. 
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The analysis also shows that the onset of the Schottky retarding region for cylindri- 
cal diodes does not coincide, in general, with the onset of monotonically retarding 
potentials. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 9, 1968, 

129-02-01-05-22. 
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APPENDIX A 


SYMBOLS 


B o 

space -charge parameter 

U r 

thermal velocity (eq. (B3)) 


(eqs. (Bll) and (B12)) 

V 

potential energy 

E 

total energy in r, 9 plane (eq. (5)} 

V 

aximuthal velocity component 

I 

electron current (eq. (1)) 


(eq. (B3)) 

I 0 

emitted (or saturation) electron 

w 

axial velocity component (eq. (B3)) 


current 

z 

axial coordinate (eq. (Bl)) 

k 

Boltzmann constant (eq. (3)) 

e o 

vacuum permittivity (mks units) 

L 

angular momentum (eq. (4)) 


(eq. (B12)) 

l 

length of diode (eq. (B5)) 

V 

parameter (eq. (B3)) 

m 

electron mass (eq. (4)) 

6 

aximuthal coordinate (eq. (Bl)) 

n 

electron density (eq. (B6)) 

£ 

transformed coordinate (eq. (B13)) 

P 

parameter, space -charge 

<P 

dimensionless potential (eqs. (1) 


integral (eq. (C30)) 


to (3)) 

P 

parameter, space -charge 

<P C 

collector potential 


integral (eq. (B9)) 

^cm 

collector potential for potential 

Q 

parameter, space -charge 

with zero slope at collector 


integral (eq. (C30)) 

^sat 

collector potential for potential 

q 

parameter, space -charge 

with zero slope at emitter 


integral (eq. (BIO)) 

^sch 

collector potential for first 

r 

radial coordinate (eq. (4)) 

potential in SRR 

S 

space -charge integral 

Subscripts: 


(eq- (Cl)) 

e 

envelope, Wheatcroft plots 

T 

emitter temperature (eq. (B3)) 

m 

potential minimum 

u 

radial velocity component 

max 

maximum value 


(eq. (B3)) 

min 

minimum value 

u . u 

upper and lower limits of 



o — o 

integration over u Q (eq. (B7)) 

Superscript: 



A 

dimensional variable 
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APPENDIX B 


MATHEMATICAL FORMULATION 

Except for notation, this appendix reiterates appendix 1 in Wheatcroft’s paper (ref. 2). 
All symbols are defined in appendix A. 


Motion of Electrons in Cylindrical Geometry 

The coordinate system to be employed is depicted in figure 9. The velocity compo- 
nents are given by 


u = 


dr 

dt 




v = 


r dg 
dt 




w = 


dz 

dt 


(Bl) 




In this system, the dimensionless radial velocity component u(r) satisfies the equa- 


tion 


u 2 (r) = <p( r) + u 2 + t?v 2 


(B2) 



Figure 9. - Cylindrical 
coordinate system. 
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where 



v = 


v 



r = 


r 



>- 


(B3) 


r , -1 „-2 

T) = 1 - r 


u 


T 



and the subscript o refers to conditions at r = r Q (r = 1). 

With this notation, the Maxwellian velocity distribution function of the emitted cur- 
rent I can be written as 

21 - 

dl = — - u e V 0 0 0 'du dv dw (B4) 

o o o o 

7 r 

A A 

where, as in equations (B3), w = w/u T> 

Electron Space Charge 

The contribution to the dimensionless electron density n(r) = n(r)/n Q of dl reach- 
ing r is given by (after integrating over w from -°° to °°) 

dn(r) = — — du dv (B5) 

t rl ru(r) ° 


where n Q is the density of emitted electrons given by 
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n 


o 


I o VtT 

2n Q l eu r 


and u(r) is given by equation (Bl). 

Integrating over u Q , v Q results in 



— o 


(B6) 


(B7) 


It should be noted that contributions to the density n(r) at any radius r will result 
from electrons flowing to the collector as well as electrons returning to the emitter. 
These various contributions are implicit in the limits of equation (B7). Their unique 
formulations are presented in appendixes C and D where u Q and u q are explicitly 
defined in terms of v . 

The integration over du Q in equation (B7) may be performed in closed form by first 
transforming the variable of integration to u(r) (eq. (B2)): 


n(r) = 


OO 


fp( r) / -r~ 2 v 2 / \ 


: / e °(erf p - erf qldv 

(B8) 

rVr/ ^ ' 

•'-co 


P 2 = u 2 + W 2 0 + <p(r) 

(B9) 

q 2 = u^+ T)v 2 + ^(r) 

(BIO) 

Poisson's Equation 



Poisson’s equation in dimensionless form for the cylindrical case becomes 

= B n(r) (Bll) 

r dr \ dr/ ° 


where B Q (in mks units) is 
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eVo 

B = 0 0 

0 e kT 


(B12) 


and e Q is the permittivity of free space. Employing the transformation, 

£ = In r 

Poisson's equation (Bll) can then be expressed as 

<P"(0 = e 2 ^ n 


(B13) 


(B14) 


where the primes denote differentiation with respect to In terms of equations (B8), 
(B9), and (BIO) become 


n(£, 



e” z (erf p - erf q) dz 


(B15) 


where 


and 


z3V o e 


P 2 = U^+ (e 2 ^ - l)z 2 + cpti) 
q 2 = u^+ (e 2 ^ - l)z 2 + 


(B16) 


(B17) 

(B18) 
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APPENDIX C 


RETARDING FIELD 
Space -Charge Integral 

Although the numerical solutions will ultimately be performed by employing the in- 
dependent variable ij and equations (B14) to (B18), in this section the independent vari- 
able r (eqs. (B8) to (BIO)) will be used purely for notational convenience and clarity 
because it precludes the need for writing the troublesome and nonpysical e^ for r. 
Furthermore, it will be convenient to define a factor S such that 


n(r) = 


2e < ? (r) 

rVff~ 


S(r) 


(Cl) 


or 


S = 



2 , 2 
-v o / r 


(erf p - erf q)dv Q 


and from equation (BIO), 


P 2 = % + -r]v 2 0 + cp( r) 
q 2 = t]Vq+ cp(r) 


(C2) 


(C3) 



Figure 10. - Typical potential. 
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A typical dimensionless potential distribution is shown in figure 10 (note by definition of 
cp , V max corresponds to The retarding region of this potential is specified by 

r < r m - It is now necessary to determine the limits "u 2 and u^. They will, in general, 
be functions of v . 


Case 1: Wheatcroft Case 

2 2 

The relation between u Q and v Q is most easily seen by a Wheatcroft plot. The 
actual plot employed by Wheatcroft (ref. 2, fig. 7) is shown in figure 11 with the notation 



Figure 11. - Wheatcroft plot. 


of the present report. Before discussing the interpretation of the Wheatcroft plot, it will 
be helpful to describe its derivation. From equation (B2), the radial velocity at any r 
of an electron emitted with initial velocity components u Q and v Q is given by 

u 2 (r) = u 2 + ?7V 2 + cp( r) (C4) 

2 2 

If a given electron (i.e. , a given ordered pair (v*,u ~)) reaches a given radius r, its 

O O n 

radial velocity there must be greater than or equal to zero. Setting u^(r) = 0 in equa- 

2 9 

tion (C4) results in a linear relation between u^ and v£: 

u 2 = -<p( r) - r?(r)v 2 (C5) 
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2 2 

The locus of points (v„,u ) satisfying equation (C5) is drawn in figure 11. None of the 

O O 9 9 

electrons whose ordered pairs (v^, u^) lie below this line can reach radius r. The con- 
verse is not true, as shown in the next paragraph. 

In figure 12 equation (C5) has been plotted for two radii. The line labeled r m cor- 
responds to the location of <p . As described previously, none of those electrons with 

2 2 2 2 
ordered pairs (v*, u Q ) in the region defined by = 0, v^ = 0, and 

2 2 

Uq < - <p(r m ) - t J( r m ) v o reach r m> Neither will all^ the electrons lying outside this 
region reach r . In fact, the shaded area in the figure shows a group of electrons out- 



side this region which do not reach r x < r m - This is portrayed in a different manner by 
the EODP given in figure 3(a). There the critical potential energy was not V max but 

; generally V occurs at a smaller radius than The correlation between 

figure 3(a) and figure 12 can best be understood by noting that a given vz in figure 12 
corresponds to a particular in figure 3(a). Increasing u* in figure 12 corresponds 
to an increase in E in figure 3(a). 

The curve in figure 12 is simply the envelope of the curves in equation (C5) for 

2 2 

r < r < r . All electrons with (v , u ) outside this envelope reach the collector, or, 
o m o o ~ 

equivalently, all electrons with (V^,E) so that E > V max in figure 3(a) will reach the 
collector. 

2 2 

The equation for the envelope (Vg,Ug) = 0 can be derived in parametric form from 
equation (C5) by the usual method to give 
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and 


v»( t) . - 

e 2 


Ug(t) = -<p( t) - 7?(t)Vg 


1< t < r 


m 


1< t < r 


m 


(C7a) 

(C7b) 

(C8a) 

(C8b) 


(the Wheatcroft case (ref. 2) explicitly requires the existence of a potential minimum in 
the interelectrode space). 

o 

From equations (C2) and (C3), the minimum radial velocity u Q at the emitter 
such that an electron will reach radius r is (see fig. 11) 


fe’“e) t -r = <°> ‘ O 

' /L x m 


u^ = -<p( r) - 7?(r)vJ 


2 2 / 2 \ 
u^ = u (v 1 
— o e \ o) 


0 < v < v (r) 
o e 


v 0 (r) < v Q < 00 


(C9) 


The maximum radial velocity for electrons passing radius r once is 


u o = °° 


0 ^ v Q — 00 


and for electrons passing radius r twice is 

= u e ( v o) ° s v„sv„(r) 


o e' 


(ClOa) 


(ClOb) 


Substituting equations (C9) and (CIO) into equation (C3) results in, for electrons passing 
radius r once, 


erf p = 1 


and for electrons passing radius r twice, 


0 < v < 00 

0 


(Cl la) 


erf p = erf p 


0 - v o - v e (r) 


(Cl lb) 
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erf q = 0 


(C 12a) 


o S v 0 s v e (r) 


erf q = erf p 


v g (r) < v Q < 00 


where 


Pe =u ef?) +T?(r)v o + ^ (r) 

Substituting equations (C 1 1) to (C13) into equation (C2) results in 


f" -v 2 /r 2 f 
! X<r) = / « ° <*<, + / 

Jq Jo 


^ -V^ 

e erf dv - 

e o 


/>« 

\(r) 


(C 12b) 


(Cl 3) 


2 / 2 

V r 

e erf p e dv Q 


r s r 


m 


or 


S,(r) = r 


£*/ 


V » <r) -vV 

e 0 erf p dv 


/ OO 

v e (D 


2 / 2 

-v Jr 

e erf p e dv Q 


r :£ r 


m 


(C14) 


where the subscript 1 refers to case 1, and v^(r) is, from equation (C7a), 

v 2(r) = -£V&± 


(Cl 5) 


Case 2: Modified Wheatcroft Case 

Wheatcroft (ref. 2) considered only the situation where r m < r c (radius of the col- 
lector). The situation that exists when the potential is monotonically decreasing with 
negative slope at the collector is a simple extension which can be completely described 
by reference to figure 13. Here the envelope jvg(t), u^(t)j = 0 proceeds in the same way 
as in case 1 from t = 1 to t = r c - Then, however, it deviates from the previous case 
(dashed curve in fig. 13) and proceeds along the straight line defined by 


u Q = -<?(r c ) -??(r c )v 0 
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Figure 13. - Modified Wheatcroft plot. 


Hence, the only change in analysis for case 1 is an obvious modification to the definition 
of the envelope, equation (C7). 


Case 3: Effective Accelerating Field for r 0 <r<r m 

Case 3 can be defined by either figure 14 or the EODP of figure 4(b). Figure 3(a) 
shows that case 1 tacitly assumes a net retarding field at the emitter. This is most 
easily observed by noting the slope of V at r Q ; a positive slope implies a retarding 
field and vice versa. Figure 4(b), however, shows case 3 to have an accelerating field in 
front of the emitter. 

2 2 

The only physically meaningful region in figure 14 is that for which v Q > 0, u Q > 0. 
The point (v^, u Q |= rvpOj indicated in this figure corresponds to r^ in figure 4(b); 
that is, an electron emitted with E = 0 and L = r Q Vj will reach radius r ^ before being 
reflected by V to the emitter. 

Figure 14 implies a collector potential cp„ well within the SRR, for it is readily 

c / 2 2\ 

apparent that any electron with emitted velocities such that (v Q , u Q I lies on or above the 
line labeled r c will reach the collector. 

For this case, where r = rj 

u^ = -cp( r x ) - 0 - v 0 - V 1 


“o = ° 


v l- v o 
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Initial tangential velocity squared, vjj 

Figure 14. -Wheatcroft plot: effective accelerating region. 

and for electrons passing rj once, 


“o = °° 

8 

VI 

o 

> 

VI 

o 

and for electrons passing r j twice, 


2 

u 0 = - <p( r c ) - V(r c ) 

o 

> 

VI 

o 

> 

VI 

o 

CM O 
> 

Hence, for electrons passing r^ once 


erf p = 1 

8 

VI 

o 

> 

VI 

o 

for electrons passing r^ twice 


erf p = erf p c 

0 ^ v ^ v 
o c 

erf q = 0 

0 - v 0 - V 1 

erf q = erf q 

8 

VI 

c 

> 

VI 

> 




r 

J 0 


C " v o /r l 
e 0 1 erf p r - 


f 


2 / 2 
-v /r 

e erf q 


(Cl 8) 


p c s - J? ( r c ) ] v o + ^ r l) " 

4-^ 0 + 


(C19) 


where from equation (C5) for = 0, 


<p( r 


v = - 


c) 


V, = - 




f?(ri) 


(C20) 


General Case for Monotonically Retarding Potentials 

Between cases 1, 2, and 3 there are a number of cases which are merely variations 
on the same theme. Rather than describe all possible cases individually, the equations of 
a general situation will be defined. It will then be shown how these equations may be 
easily specialized to represent all possible cases. 

The EODP for the general case is shown in figure 3(d). The Wheatcroft plot is 
shown in figure 15. In this figure, the envelope of interest is given by the implicit equa- 
tions 

jv^t), U 2(t)J = 0 l<t<r a (C21a) 

JV|(t), uf(t)J = 0 r b -t-r m (C21b) 

p 

Let v~, be defined by the point of tangency between the envelope and equation (C5). Then 
r 2 2 

for r such that v < v_, figure 15 can be represented by figure 16 where the envelope 
/ 2 2\ x r 

( v e’ u e) =0 is S ivent) y 
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( v e- u eH v l' u l) =0 

M-M-o 


V > V 
O X 


0 < v < v 

O X 


(C22) 


Now figure 16, with the envelope defined by equation (C22), represents the same situation 
as figure 11, case 1. Hence, from case 1, 

r /% (r) - V 2 / r 2 C °° -v 2 /r 2 

S(r) = ^ + j e 0 erf p g dv Q - j e 0 erf p g dv Q v g (r) > v x 

Jq Jy ( r ) 


L 


°° 2/2 

-v_/r 

e erf p g dv Q 


v e (r)>v x 


(C23) 


and p e is given by equation (Cl 3). 

The next region of interest (see fig. 15) is that for which v g (r) = v^(r) < v x and 



Initial tangential velocity squared, 
Figure 17. - General case: regime 2. 


r < r & . The Wheatcroft plot for this case is shown in figure 17. Hence, 


u^ = -cp( r) - 7j(r)Vg 


0 < v < v 
o r 


2 . 2 / 2 ' 


v — v 
o r 
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for electrons passing r once, 


U 2 = oo 


for electrons passing r twice, 


—2 2 / 
u o = u 2 (v o } 


0 < y — oc 

o 


0 - v o - v x 


and for electrons passing r once, 


erf p = 1 


0 < V <°o 

o 


(C24a) 


for electrons passing r twice 


erf p = erf p 0 < v < v 

*e ox 


erf q = 0 
erf q = erf qj 


0 < v < v 
o r 


V 2: v 
o r 


(C24b) 


(C25) 


where p is given by equation (Cl 3) since u 9 (r) = utr) of equation (C22) in the range 
0 ^ v Q /v x , and 

qj = Uj(r) + T7(r)v^ + <p(r) (C26) 

The space -charge integral for this region then becomes 


S(r) = 


U f ' X e - V ° /r2 
4) 


/ 


X -v^/r 2 


erf p e dv o " / e erf q l dv o 


/ 


2 / 2 
-v /r 

e erf q x dv 0 v r < v x , r < r & (C27) 
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Regime 3 is shown in figure 18; for this regime 

2 2 
= -<p( r) - 7?(r)v 0 


0 < v_ < v 

o r 


% = "fa 


V > V 

o r 


for electrons passing r once, 


U 2 = CO 


0 < Vq ^ 00 


for electrons passing r twice, 


*o - 


for electrons passing r once, 


erf p = 1 


0 < v < v 
o x 


0 < v ^ °° 

o 


for electrons passing r twice, 


erf p = erf p 


0 s v o - v x 


erf q = 0 


0 < v Q C v i 
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erf q = erf 


q l v o > v r 


and 


S(r) = + 

2 


jf 


X " v o /r2 
e erf p e dv 


.-jp.-s- 


erf q l dv o 


/ 


2 / 2 
-v o /r 


erf qj dv c 


v r < v x , r > r a (C28) 


where p g is given by equation (C13) and qj by equation (C26). 



p 

Initial tangential velocity squared, v£ 
Figure 19. - General case: regime 4. 


The final regime of interest in the general case is shown in figure 19. The analysis 
of this regime is identical to that of regime 1 (fig. 16). Hence, from equations (C13) 
and (C20) to (C23) 

2 / 2 
-v_/r 

e erf p e dv Q v e ( r ) < v x ( C29 ) 
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The space -charge integral for the general case, equations (C13) and (C20) to (C28) 
may be put in the following more concise form : 

•B z*°o /%B 


S(r) = 


f - 0 .r Pd v o - r 

Jq Jo Ja 


Q dv 


2 / 2 

2 “ v o/ r 
= e erf 

2 / 2 

2 _v n / r 
= e 0 erf 


u^(r) + 7?(r)v^ +^(r) 


Uj(r) + rj(r)v^ + <p(r) 


v%) = - 
1 2 

Uf(t) = - <?(t) - 7?(t)v? (t) 

where i = 1 for 1=£ t < r , and i = 2 for r b ^ t ^ r m> 




v l >v x 


•A 


v 2 <v x 


(C30) 


(C31) 


(C32) 


(C33) 


>j)-M 

and the limits A and B are defined in table I. Note that from (C32), = - jr^’(r)J/2. 

The Wheatcroft case (case 1) corresponds to v = 0. Case 3 is obtained for r = r 

with u., = 0, 0 < v < oo. 

1 ’ o 


TABLE I. - LIMITS A AND B OF SPACE - 
CHARGE INTEGRAL 


Regime 

Criteria 

A 

B 

1 

V > V 

V 

V 


r x 

r 

r 

2 

v < v , r < r 

V 

a 

V, 


r x a 

r 

X 

3 

V r < V r > r a 

bv int (r) 

v x 

4 

V < V 


v 


r x 

v r 

r 


intersection point between v, (t) and v 0 (t). 
b i z 

Intersection point between 

u o = ' <p(r) - J?(r)v^ and Uj = Uj(r). 
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General Case for r c <r m 


After modifying v 2 (t) (eq. (C32)) as indicated in case 2, equations (C30) then 
(C33) with table I represent the most general situation. 


Current to Collector 

To obtain the current to the collector, it is only necessary to integrate equation (B4) 
along the envelope , u^j of the generalized Wheatcroft plot figure 15 (after integrating 
over dw n ): ' 

2 2 
- v" 

e e °dv Q (C34) 
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APPENDIX D 


ACCELERATING FIELD 


An accelerating field is present for r > r m (see fig. 10). The Wheatcroft plot for 
the general case (r < r c ) is shown in figure 20. 

The analysis is the same as that for a retarding field with 


2 2/2 
Ho - u e\o 


— 2 = 00 
u o 


where = u g^ v e) is defined b y equation (C32) in the general case. Hence, 


/ °° 2/2 . 

e 0 erf qK) dv 0 


where 


q = u e( v oj + 7 ? ^ v o + ^ 


v u o - -<P m - 


-u l- 


Initial tangential velocity squared, 

Figure 20. - Wheatcroft plot for accelerating region. 
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